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ABSTRACT 

We present deep multiband (F435W, F625W, and F658N) photometric data of the Globular Cluster 
uj Cen collected with the Advanced Camera for Surveys on board of the Hubble Space Telescope. 
We identified in the (F435W-F625W, F435W) plane more than two thousand White Dwarf (WD) 
candidates using three out of nine available pointings. Such a large sample appears in agreement with 
predictions based on the ratio between WD and Horizontal Branch (HB) evolutionary lifetimes. We 
also detected « 1600 WDs in the (F658N-F625W, F625W) plane, supporting the evidence that a large 
fraction of current cluster WDs are Ha bright. 

Subject headings: globular clusters: general — globular clusters: omega Centauri 



1. INTRODUCTION 

Among Galactic Globular Clusters (GGCs), lu Cen is 
the most massive one (M = 5 x 10^ Mq, Meylan et 
al. 1995) and the only one which clearly shows a well- 
defined spread in metallicity. According to recent es- 
timates based on sizable samples of evolved Red Giant 
Branch (RGB) and sub-giant branch stars, the metallic- 
ity distribution shows three peaks around [Fe/H] = —1.7, 
— 1.5 and —1.2 together with a tail of metal-rich stars ap- 
proaching [Fe/H] « —0.5 (Norris et al. 1996; Hilker et al. 
2004). During the last few years it has also been sug- 
gested that UJ Cen harbors multiple stellar populations 
(Lee et al. 1999) characterized by different ages (Ferraro 
et al. 2004; Hughes et al. 2004), He abundances, and dis- 
tances (Bedin et al. 2004; Norris 2004; Freyhammer et 
al. 2005). More recently, SoUima et al. (2005) suggested 
that RGB stars present five peaks due to different stellar 
populations. Moreover, Piotto et al. (2005) found, that 
the bluer main sequence (MS) detected by Bedin et al. 
(2004) is more metal-rich than the red MS. To account 
for this feature they suggested a possible variation of the 
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He content among MS stars. 

Even though lo Cen presents properties which need to 
be properly understood, its stellar content is a gold mine 
to investigate several open problems concerning stellar 
evolution and its dependence on the metallicity. This 
applies not only to evolved stars such as RR Lyrae, hot 
HB stars, and the tip of the RGB, but also to the ex- 
pected population(s) of WDs. The search for WDs in 
GGCs has already been successful, and several WD sam- 
ples have been identified (Hansen et al. 2003; Moehler et 
al. 2004). So far, the richest sample of cluster WDs (222) 
was detected in M4 by Hansen et al. (2004). The detec- 
tion of WDs in UJ Cen dates back to Ortolani & Rosino 
(1987) who selected two dozen WD candidates on the 
basis of ground-based data, and to Elson et al. (1995), 
who detected four WD candidates using HST data. In 
this paper, we present preliminary results concerning the 
identification of more than 2,200 WDs in uj Cen on the 
basis oi B,R, Ha data collected with the Advanced Cam- 
era for Surveys (ACS) on board of HST, and available 
on the HST archive. The reader interested in thorough 
reviews concerning recent theoretical and empirical re- 
sults is referred to (Fontaine, Brassard, & Bergeron 2001; 
Koester 2002; Prada Moroni & Straniero 2002; Hansen 
& Liebert 2003; Hansen 2004, and references therein). 

2. OBSERVATIONS AND DATA REDUCTION 

Current data were collected with nine pointings of 
the ACS camera across the center of the cluster. The 
3x3 mosaic covers a field of view of « 9' x 9'. Four 
images per field have been acquired in three different 
bands, namely F435W, F625W, and F658N (hereinafter 
B, R, and Ha)- Three deep (340s) exposures were 
secured for the B and R-band, respectively, while the 
exposure time for the four Ha images was 440s each. 
The nine fields were independently reduced with the 
DAOPHOTII/ALLFRAME (Stetson 1994). An individ- 
ual PSF has been extracted for each frame by adopting, 
on average, Ri200 bright isolated stars. The individual 
catalogues were rescaled to a common geometrical sys- 
tem with DAOMATCH/DAOMASTER. The final cat- 
alogue includes approximately 1.2x10^ stars, and the 
B,B — R CMD shows a well-populated cluster MS to- 
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gather with a sizable sample of stars covering a wide 
region to the left of the MS (Bono et al. 2005). 

From this catalogue we selected all the stars bluer 
than MS stars and fainter than extreme HB stars. We 
ended up with a sample of approximately 45,000 stars 
distributed over the nine pointings. Among them we se- 
lected three pointings (2, 7, and 9) that include ~ 14, 000 
WD candidates. These stars have been identified in in- 
dividual deep B,R,Ha images and we performed once 
again the photometry using ROMAFOTwo. Individual 
stars have been interactively checked in every image, and 
the WD candidates were measured either as isolated stars 
or together with neighbor stars, where crowding affects 
the photometric accuracy. Note that a substantial frac- 
tion of the stars located close to WD candidates are truly 
MS stars, i.e. they did not belong to the original sam- 
ple of stars located on the blue side of the MS. Note 
that the 80% of the selected stars revealed to be either 
cosmic rays, or spurious identifications of faint stars lo- 
cated close to saturated stars, or too faint objects to 
be safely measured on individual images. The photo- 
metric calibration was performed in the Vega System 
( |http: / /www.stsci.edu/hst/acs/documentS| ) . 

3. RESULTS AND DISCUSSION 

In this investigation we present preliminary results based 
on three out of the nine pointings. Figure 1 shows the 
current sample of w Centauri stars in the B — R,B plane. 
Interestingly enough, the original sample splits into two 
well-defined sequences, the redder one made by MS stars 
and the other one made by 2,212 blue objects that ranges 
from B « 22, B - i? « 0, down to B w 27, B - i? w 0.8, 
covering the expected region of cluster WDs. 

Using GOODS data collected with ACS in the bands 
F435W and F606W, we estimated that the expected 
number of field galaxies with 22 < i? < 26 and - 
0.2 < B — R < 0.5^^ in the same area covered by cur- 
rent observations is w 60 (Grazian et al. 2005, priv. 
comm.). The number of field stars is also negligible, be- 
cause halo and disk stars peak around B — R = 0.7 and 
B - R = 1.8, respectively (King et al. 1990). More 
detailed estimates based on radial velocity (Suntzeff & 
Kraft 1996) or proper motion (van Leeuwen et al. 2000) 
measurements in to Cen suggest that at most two dozen 
of field stars might be located inside the area covered 
by current data. Finally, theoretical Galactic models 
(Castellani et al. 2002) suggest that in the same area 
covered by current observations are present at most six- 
teen field WDs brighter than B = 25.5 and fifty-five field 
WDs brighter than B = 28. 

This evidence suggests that we are facing with a bona 
fide sample of cluster WDs, including more than 2,000 
objects, thus the largest sample of WDs ever observed 
in a stellar cluster. Data plotted in Fig. 1 clearly show 
that, thanks to the sizable sample of w Cen stars, the 
cooling sequence show up at B « 21. This bright limit, 
if we assume for lo Cen an apparent distance modulus 
DMb ~ 14.21 (Thompson et al. 2001), implies that cur- 
rent data are tracing the cooling history of cluster WDs 

Note that to perform this estimate we accounted for the 
diflference in magnitude between AB and Vega systems and for 
t he differenc e between the filter F606W and the filter F625W 
http:/ /www. stsci.edu/hst/acs/documents I. 



at least below Mb ^ 7. According to predictions by Al- 
thaus, & Benvenuto (1998) for a 0.5Mo WD with a CO 
core and pure H atmosphere models by Bergeron, We- 
semael, & Beauchamp (1995)^^ the current WD sample 
provides the opportunity to investigate the WD cooling 
for luminosities ranging from ten times the solar lumi- 
nosity down to logL/L© 3.1 {B - 27, Mb ^ 12.8). 

By adopting the same theoretical predictions one can 
easily recognize that the huge number of WDs should 
not be a real surprise. As a matter of fact, the num- 
ber ratio between WD and HB stars, for not too long 
cooling times, is simply given by the ratio of the life- 
times spent during these two evolutionary phases. On 
the basis of the shallow ACS photometry in the B, R 
bands, recently presented by Freyhammer et al. (2005), 
we found in the same three fields the occurrence of 630 
HB stars. The typical evolutionary lifetime for HB stars 
is ~ 1 — 2 X 10* year, depending on their ZAHB effec- 
tive temperature (Castellani et al. 2004), while for a 
WD with a CO core and M = 0.5Mq at Mb ~ 11.3 
(logi/Lg = -2.24) the lifetime is ~ 3.5 x 10* years. 
This means that we expect to detect roughly twice as 
many WDs than HB stars brighter than B w 25.5. Inter- 
estingly enough, we detected approximately 1200 WDs 
brighter than this magnitude limit. However, the current 
estimate should be considered as a robust lower limit, 
since preliminary artificial star experiments suggest that 
at this limiting magnitude the completeness is « 85%. 

An unexpected observed feature in Fig. 1 , is the steady 
increase in color dispersion when moving toward fainter 
WD magnitudes. In order to assess whether this spread 
in color might be due to photometric errors, we per- 
formed an empirical test. We estimated the ridge line 
of WDs and the distance in color of individual objects 
from the ridge line. Figure 2 shows the color distribu- 
tion of WDs in the magnitude interval B — 25 ± 0.5 
together with its gaussian fit (solid line). Then we esti- 
mated standard deviation in color of the same selected 
WDs (intrinsic errors) and we found ctb-b. = 0.087. The 
dashed line plotted in Fig. 2 shows the expected color 
distribution for the same sample of WDs (712) according 
to the assumption that their B, R magnitudes would only 
be affected by gaussian intrinsic errors. Data plotted in 
Fig. 2 indicate that the sigma of the gaussian fit to ob- 
served WDs is a factor of two larger than expected. On 
this basis, we estimated that the two distributions dif- 
fer at 99% confidence level. This finding, taken at face 
value, indicates that the color dispersion might be real. 

In order to investigate the WD location in the CM dia- 
gram we adopted the WD cooling sequences for CO core 
and H envelopes constructed by Althaus, & Benvenuto 
(1998). Theoretical predictions have been transformed 
into the observational plane by using pure H atmosphere 
models constructed by Bergeron et al. (1995). Pre- 
dicted cooling sequences for M = 0.5, 0.7, 0.9Mq have 
been plotted by adopting canonical estimates for cluster 
reddening {E{B - y) = 0.11 ± 0.02, Lub 2002) and dis- 
tance modulus {fj, — 13.7 ±0.11, Thompson et al. 2001). 
The reddening in the B, R, Ha bands was estimated us- 
ing the extinction model of Cardelli et al. (1989). 

The comparison between theory and observations (see 
panel a) of Fig. 3) discloses a good agreement for the 

See also http: //www. ASTRO.UMontreal.CA~bergeron/CoohngModels/| 
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bright portion (22 < B < 23). However, the theoretical 
sequences toward fainter magnitudes appear to fit the 
blue (hot) edge of the observed WDs. It is noteworthy, 
that a WD with a stellar mass M = 0.5Mq is a lower 
limit for WDs with CO-core, actual WDs are expected 
to be slightly more massive (~ O.SSA/q, Renzini et al. 
1996). Current uncertainties on reddening (0.02 across 
the cluster, Schlegel et al. 1998) and on cluster distance 
cannot account for the observed systematic drift in color. 
This evidence suggests that the observed WD cooling 
sequence is redder than expected. 

In principle one can find several plausible reasons for 
such an occurrence. In particular, we note that for 
22 < i? we are already below the so-called "DB gap" 
(Hansen & Liebert 2003), i.e. WDs with He atmospheres 
could be present. Data plotted in panel b) of Fig. 3, show 
that by adopting WD coohng sequences for CO core, and 
He envelopes by Benvenuto, & Althaus (1997) together 
with He atmosphere models by Bergeron et al. (1995), 
the cooling sequence is indeed moving toward redder col- 
ors. Note that the cut-off of the cooling sequences in the 
bright region is due to the fact that He atmosphere mod- 
els do not cover this temperature region. The comparison 
between theory and observations indicates that the bulk 
of the observed WDs might be of the DB type. This 
evidence, as suggested by the referee, is in contrast with 
current spectroscopic findings by Moehler et al. (2004) 
concerning cluster WDs. They identified of H Balmer 
lines in four and seven WDs in NGC 6397 and in NGC 
6752, thus suggesting that these objects are of the DA 
type. The same outcome applies to field WDs for which 
DA WDs are « 80% of the entire sample (Koester & 
Chanmugam 1990). An explanation of this discrepancy 
should await for new spectroscopic measurements. 

As shown in the bottom panel of Fig. 3, a different 
possibility is given by the occurrence of He-core WDs. 
Interestingly enough, predicted WD cooling sequences 
for He core structures by Serenelli et al. (2002), trans- 
formed into the observational plane by adopting H at- 
mosphere models, also account for the observed distri- 
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bution. This working hypothesis, once confirmed, indi- 
cates that a substantial fraction of WDs in u! Cen might 
be the evolutionary aftermath of binary systems. Ac- 
cording to current knowledge, candidate cluster He core 
WDs have only been identified in systems with high 
central densities and short two body relaxation times 
(NGC6397, logpo = 5.68 Lqpc-^, Taylor et al. 2001; 
47 Tuc, logpo = 4.77 Lqpc'^, Edmonds et al. 2001). 
However, lo Cen presents a relatively low central density 
(logpo = 3.12 Lqpc~^), and therefore fills an empty re- 
gion of the parameter space predicted by Hansen et al. 
(2003, see their figure 8). Current findings are marginally 
affected by the adopted theoretical predictions. Differ- 
ent sets of WD cooling sequences (Fontaine et al. 2001; 
Prada Moroni & Straniero 2002), transformed into the 
B — R, B plane by adopting the same atmosphere mod- 
els, agree quite well with each other. 

Finally, it is worth noting that current — -R, i? data 
(see Fig. 4) , show that a good fraction of detected WDs 
1,600) are Ha bright, and indeed they attain Ha — R 
colors systematically bluer than predicted. We would 
like to mention that blending and/or binarity could also 
affect observed colors. Therefore, no firm conclusion can 
be reached on the basis of current data. More detailed 
photometric and spectroscopic investigations of this large 
sample of cluster WDs will provide fundamental hints 
concerning their cooling and progenitors. 
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Fig. 1. — Color-Magnitude diagram in the B,B — R bands. The number of WD candidates are labeled. Error bars display intrinsic 
photometric errors. 
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Fig. 2. — Color distribution of WDs in the magnitude interval B = 25 ± 0.5. The solid line shows the gaussian fit to the observed 
distribution, while the dashed line the expected distribution for the same sample of WDs in the case that their colors would only be affected 
by gaussian photometric intrinsic errors. 
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Fig. 3. — Color-Magnitude diagram in the B, B — R bands. Panel a) - Solid, dotted, and dashed lines show predicted WD cooling 
sequences for CO core and H envelopes with stellar masses equal to 0.5, 0.7, and 0.9 M/Mq (Althaus, & Bcnvenuto 1998). Theoretical 
predictions have been transformed into the observational plane by adopting H atmosphere models. Diamonds show a 12 Gyr isoclirone for 
Z=0.001 and Y=0.232 constructed by Cariulo et al. (2004, M > 0.5Mq) and by Baraffe et al. (1997, M < 0.5Mq). Panel b) - Same as 
panel a), but for WD cooling sequences, with CO core and He envelopes (Benvenuto, & Althaus 1997), transformed into the observational 
plane by adopting He atmosphere models. Panel c) - Same as panel a), but for WD cooling sequences with He core and H envelopes 
(Serenelli et al. 2002), stellar masses equal to 0.30, 0.39, and 0.45 M/Mq and H atmospheres. 
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Fig. 4. — Color-Magnitude diagram in the R, Ha — R bands. Predicted cooling sequences are the same as in Fig. 4. The thick dashed 
line marks the detection limit. Error bars display intrinsic photometric errors. 



